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Abstract. We present an algorithm to compute a primary decomposition 
of an ideal in a polynomial ring over the integers. For this purpose we use 
algorithms for primary decomposition in polynomial rings over the rationals 
resp. over finite fields, and the idea of Shimoyama-Yokoyama resp. Eisenbud- 
Hunecke-Vasconcelos to extract primary ideals from pseudo-primary ideals. A 
parallelized version of the algorithm is implemented in Singular. Examples 
and timings are given at the end of the article. 



1. Introduction 

Algorithms for primary decomposition in Z[xi, . . . , x n ] have been developed by 
Seidenberg (cf. [Se]), Gianni, Trager, Zacharias (cf. jGTZj ) and Ayoub (cf. [A]). 
Within this article we present a slightly different approach which mainly uses pri- 
mary decomposition in polynomial rings over a field and therefore seems to be 
much more efficient. In particular, it uses primary decomposition in Q[aci, . . . , x n ] 
resp. F p [xi, . . . , x n ] as well as the computation of the minimal associated primes of 
an ideal in ¥ p [x\,... , pseudo-primary decompositiorfl, and the extraction of 
the primary components. The essential difference compared to the corresponding 
algorithm proposed in [GTZ] is as follows: the primary decomposition of an ideal 
I in Z[zi, . . . , x n ] with I n Z = (q) such that q ^ is obtained by computing the 
minimal associated prime ideals of I¥ p [x\, . . . ,x n ] for all primes p dividing q and 
extracting subsequently the primary ideals. 

Let x — . . . ,x n } always denote a set of indeterminates and let I C Z[x] 
be an ideal. We use the following known facts from commutative algebra for our 
algorithm: 

(1) If I n Z = (0), then there exists an/ieZ such that J : h = IQ[x] H Z[x] 
and I = (I : h) H (I, h) (cf. [Si], Theorem 2). 

(2) If / n Z = (0) and IQ[x] = Q 1 D . . .DQ S is an irredundant primary decom- 
position with Pi = Jq, then iq[x] nz[i] = (Q x nz[i]) n . . . n (Q s n z[x}) 



is an irredundant primary decomposition and Pi flZ[x] = y Q i n Z[x] (cf. 
[Se] . Theorem 3). 
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1 One can choose one of the modern algorithms, cf. |DGPI . IEHVI . IGTZI . \§Y] . 

2 An ideal is called pseudo-primary if its radical is prime, cf. [EH V) . |SY| . 
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(3) If / n Z = (q) such that q 7^ and q = p^ 1 ■ ■ -p v r T with pi, . . . ,p r pairwise 
different primes, then / = C\l = i(I,Pi % ). 

(4) If If] Z = for some prime p and Pi, . . . , P. are the minimal associated 
primes of 7F p [x], then the canonical liftingsj Pi,...,P s to Z[x] are the 
minimal associated primes of I. 



If v = 1 let 7Fp[x] = <3i H . . . H Q s be an irredundant primary decompo- 
sition with associated primes Pi, . . . , P s and Q%, . . . , Q s , Pi, . . . , P s be the 
canonical liftings to 1[x] . Then / — Q\ n . . . n Q s is an irredundant primary 
decomposition with associated primes P%,...,P S . 
The following result can easily be adapted to Z[x]. 

(5) If P is a minimal associated prime of /, then / + P m is a pseudo-primary 
component of / for a suitable m 6 N, i.e. the equidimensional part of 
j _l_ pm j g p r j mar y component of J associated to P. For any m let Q m 
be the equidimensional part of I + P m . Q m is a primary component of / 
with associated prime P if Q m = IZ[x]p fl Zfx] (cf. [EHV] ) - 

Alternatively we can compute a separatoiQ s of I w.r.t. P and obtain 
by I : s°° a pseudo-primary component of I (cf. |SY] ) . 

(6) If Q\, . . . ,Q S are the primary components of I associated to the minimal 
associated prime ideals and J = Q± D . . . (1 Q s , then there exists a natural 
number m such that I = J n (7 + (I : J) m ). 

Consequently, by applying (l)-(6), we can reduce the computation of the primary 
decomposition in Z[x] to the computation of the primary decomposition in Q[x], the 
computation of the minimal associated primes in ¥ p [x] , and the extraction of the 
primary components in 7L\x\. In this connection, the extraction has to be generalized 
to polynomial rings over principal ideal domains (cf. Lemma l2.5[) . In section [2] we 
state the results used in the algorithm, whereupon in section [3] we explain our 
algorithm which has been implemented in Singular in a parallel version. Finally 
we give some examples and the corresponding timings in section |4j 



Definition 2.1. Let I C 1[x] be an ideal and > be a monomial ordering on Z[x]. 
A subset G C I is called a Grobner basis of I w.r.t. > if the leading ideal of G 
equals the leading ideal of I. G is called a strong Grobner basis if for all / € I 
there exists ag&G such that LT( 5 )| LT(/)0 

Lemma 2.2. Let G = {gi, . . . , gt} Q Z[x] and I = (G)Z[x]. Assume that If] 
Z = (0) and G is a Grobner basis of I<Q[x] w.r.t. some ordering. Let h — 
lcm(LC((7i), . . . , LC(gfc)) be the least common multiple of the leading coefficients 
of g u ...,g k - Then IQ[x] n Z[x] = I : h°° . Moreover, if I : h°° = I : h m for some 
natural number m, then I = (I : h m ) n (I, h m ) . 

The proof of Lemma 12.21 is similar to the corresponding proof for polynomial 
rings over a field (cf. |GPj . Proposition 4.3.1). 

■^Choose generators in V p [x] and lift the coefficients to non-negative integers smaller than p. 
^We call s a separator of / w.r.t. P if s £ P and s is contained in all other minimal associated 
primes of J. 

^We use the notations of |GP| for the basics of Grobner bases. Especially LT(/) denotes the 
leading term (leading monomial with leading coefficient) of / w.r.t. the ordering >. The theory 
of (strong) Grobner bases over principal ideal domains can be found in IAL| . section 4.5. 



2. Basic definitions and results 
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Remark 2.3. The saturation / : h°° can be computed in Z[x] similarly to the case of 
a polynomial ring over a field by computing a Grobner basis of (I, Th — 1) Z[x, T] 
w.r.t. an elimination ordering for T: 

L : h°° = (I, Th-l) Z[x, T] nZ[i], 

A natural number m satisfying / : h°° = I : h m can be found by computing the 
normal form of h l g w.r.t. I for each generator g of I : h°° and increasing I € N. 
More precisely, if the normal form of h l g w.r.t. / is zero for each generator g of 
L : h°° then h l ■ (I : h°°) C I, i.e. I : h°° = I : h l . 

Lemma 2.4 (cf. SY J). Let I C Z[x] be an ideal with more than one minimal 
associated prime, P a minimal associated prime and s P a separator, i.e. s is 
contained in all other minimal associated primes of I. Then I : s°° is a pseudo- 
primary component of I , and s can be chosen as 

n 

QGminAss(J) 

where sq is an element of a Grobner basis of Q which is not in P. 

Lemma 2.5 (Extraction Lemma, cf. [GTZj ). Let L — Q n J be pseudo-primary 
with yfl = P and Q be P -primary with ht(Q) < ht( J). Let P n Z = (p) for some 
prime p and u C x be a maximal independent set of variables for P = PF p [x@. Let 
R := ZfttbjA, then the following hold: 

(1) LR[x \ it] H Z[x] = Q 

(2) Let G be a strong Grobner basis of L w.r.t. a block ordering satisfying 
x \ U 3> u. Then G is a strong Grobner basis of IR[x \ it] w.r.t. the 
induced ordering for the variables x \ u. 

(3) Let G = {gi,...,g k } be as in (2), LT R [ XX „] (g;) = p u 'a t (x \ it) ft with 
ai G Z[u] \ (p) for i = 1, . . . , k, and h — lcm(ai, . . . , a^). Then LR[x \ it] f~l 
Z[x] = L :h°°. 

Proof. 

(1) Let K = y/J imdK = K¥ p [x] then ~K D P = P¥ p . This implies that 
if n F p [it] 7^ (0) since u C a; is maximally independent for P and therefore 
K H (Z[it] n (p)) ^ 0. Thus it holds JP[a; \ it] = \ u]. Finally, because 
Q is primary, we obtain LR[x \ u] n Z[x] = QR[x \ u] n Z[x] = Q. 

(2) Let / G iP[x n u] and choose s G Z[u] \ (p) such that s/ G I. Since G 
is a strong Grobner basis of I there exists a g £ G such that LT z r,,,i(<7) | 
LT z r a .i(s/). As a polynomial in i \ w with coefficients in R, the element 
s/ can be written as sf = p v a(x \ it) Q + (terms in x \ it of smaller order) 
with a G Z[u] \ (p). If p T is the maximal power of p dividing the leading 
coefficient LCz[ x ] (<?) of g then r < v since LTg^j (s/) = p" LT^ (a)(a;\ii) Q . 
Now we can write 5 as an element of R[x \ u] w.r.t. the corresponding 
ordering, i.e. g = p^b(x \ u) 13 + (terms in x \ it of smaller order) with 
b G Z[u] \ (p) and p < r < t/. By definition we have LT^u^] (5) = 
p^b(x\u)P resp. LTjj[ xxu ](/) = p"^(x\u) a and on the other hand it holds 
LTz[x](5) = P AI LT z[2 .](6)(a; \ u)^ resp. LT Z [ x] (s/) = p" LT Z [ x] (a) (a; \ u) a . 

®u C i' is called a maximal independent set of variables for P C ¥ p [x] if P n F p [it] = {0} and 
#u = dim(¥ P [x]/P); cf. \UP] . 
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Thus the assumption LTgu.1 (g) | LT%[ x ](sf) implies (x \ u)P \ (x \ u) a and 
consequently LTr[ x ^ u ] (d) I LT^^] (/). This proves (2). 
(3) Follows from (2) similarly to the proof for fields (cf. [GTZ] . [GT] ). 

□ 

The following Lemma is a consequence of the Lemma of Artin-Rees (cf. |GPj ). 

Lemma 2.6. Let I C Z,[x] be an ideal and J the intersection of all primary compo- 
nents of I associated to the minimal prime ideals of I. Then there exists a natural 
number m such that I — J (~\ (I + (I : J) m ). 

Notation 2.7. Given an ideal I C 7L[x\ we can always choose a finite set of polyno- 
mials Fi = {fx, f k ] such that / = (Fi) and we denote F ( / n) := {/f , . . . , 
for m € N. 

Corollary 2.8. With the assumptions and notations of Lemma \2.6\ there exists a 
natural number m such that I = J PI (/ + (-F 1 }™'')). 

Proof. Due to Lemma [2.61 there exists an m such that / = J <~) (I + (I : J) m ). 
Now we have K Jn(/+ (^/ ( "/ ) >) C J n (I + (I : J) m ) = / and therefore 
7 = Jn(/+(F / ( ™ ) )). □ 

Remark 2.9. The corollary is very important from a computational point of view 
because {Fjjj) has fewer generators than (I : J) m . 



3. The algorithms 

In this section we present the algorithm to compute a primary decomposition of 
an ideal in a polynomial ring over the integers by applying the results of section [5] 
resp. the introduction (section [1]) . 

Algorithm Q] computes the primary decomposition of an ideal in Z[a;Q with the 
aid of algorithms [5] and [3] which we introduce subsequently in detail. 

Remark 3.1. Algorithm Q] can easily be parallelized by computing - depending on 
the prime factorization q — p 1 ^ 1 ■ ■ ■ p v r r of q where (q) = I PI Z - either the primary 
decomposition or the set of minimal associated primes in positive characteristic in 
parallel. If = 1 we have to compute the primary decomposition whereas, if i>i > 1, 
we have to compute the minimal associated primes of I¥ Pi [x] in F Pi [x] . These r 
computations in positive characteristic are independent from each other such that 
they can also run separately in parallel on at most r processors if available. 



7 The corresponding procedures are implemented in Singular in the library primdecint . lib. 
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Algorithm 1 primdecZ 



Input: Fj = {fx, . . . , /&}, I = (Fi)Z[x], optional: a test ideal T. 
Output: L := {(Qx,Pi), ■ ■ ■ , (Qs,P s )}, I = Qi H . . . fl Q s irredundant primary 
decomposition with Pi = \fQl. 

if T is not given in the input then 

T:= (1); 
G := strong Grobner basis of I; 
q := generator of 7 n 
if q = then 

compute h € Z such that I : h = IQ[x] nZ[i]| 

compute Q l} . . . ,Q S , an irredundant primary decomposition of IQ[x] and Pi = 
\j~Qi the associated primes; 

compute Qi = Q l n Z[x], Pi = Pi n Z[x]0 
L:={(Q 1 ,P 1 ) ) ...,(Q S ,P S )}; 

A/ := PRIMDECZ((7, h)) & remove redundant primary ideals from M; 
return LU M; 
else 

compute = p^ 1 . . . p" r , the prime factorization of q; 
for i = 1, . . . , r do 
if Vi = 1 then 

compute Li = {(Qj \ p\ ''),••■ , (Qi P^ ')}, the primary decomposition 
of/F p Jx]; 

:= {(Q«, pW), . . . , (Q^Pif)}, the lifting of to Z[x]0 

else 

compute Ai = {P 1 , .. . ,P S . }, the set of minimal associated primes of 
I¥ Pi [x] and independent sets of variables , . . . , u% > for P^,...,pif; 

A 4 := {Pf\ . . . , Pjf }, the lifting of A, to Z[x]; 
for j = 1, . . . , Si do 

:= EXTRACTZ^^^Pj 1 ',?!^); 

£< :={(Qi^f ),..., (Q^pf)}; 

L := Li U . . . U L r ; 

compute J, the intersection of all primary ideals in L and T; 
if J = 7 then 
return L; 

compute Fj.j such that (Fi-.j) = I . J; 
compute m such that J fl (J + (P/.j }) = 7; 

M := primdecZ(7+ (Pj™'), J) & remove redundant primary ideals from M; 
return LU M; 



5 is either or the unique element in G of degree 0. 
9 /i is a suitable power of the least common multiple of all leading coefficients of elements in G, 
cf. Lemma [221 

10 Qi resp. Pi are primary resp. prime due to (4) of the introduction (section [TJ. 

1J If J = (Fj) C ¥ p [x] then its lifting is obtained by (p, Fj) with the canonical lifting of Fj. 
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The algorithm to compute the separators is based on Lemma 12.41 



Algorithm 2 separatorsZ 

Input: B a list of prime ideals generated by a Grobner basis w.r.t. some ordering, 

not contained in each other, P € B. 
Output: Polynomial s such that s P, s 6 Q for all Q S B\{P}. 

for Q e B\{P} do 

choose sq in the Grobner basis of Q such that sq ^ P; 

return Uq £ b\{p} s Q5 



The algorithm to extract the primary component from the pseudo-primary com- 
ponent is based on the Extraction Lemma 12.51 



Algorithm 3 extractZ 

Input: I C Z[x] an ideal, B the list of minimal associated primes oi I, P E B with 
P H Z = (p) for some prime p, a C i an independent set of variables for PF p [x]. 
Output: The primary component Q of I associated to P. 

s := separatorsZ (P, B); 
1 = 1: s°°; 

compute G = {g\, . . . ,gu}, a strong Grobner basis of I w.r.t. a block ordering 
satisfying x \ u ^> u\ 

compute {a x ,.. .,a k } such that LC ZM (p> [ xv „] = p"* • a 4 with a, G Z[n] \ (p); 
compute = lcm(ai, . . . , afc), the least common multiple of ax, . . . , a^; 
return / : h°°; 



Example 3.2. Consider I = (9,3x, 3y), P = (3), u = {x,y} and P = {P} in Z[x,y}. 
Then we obtain s = 1, h = xy and thus 7 : h°° — (3). 

4. Examples and timings 

In this section we provide examples on which we time the algorithm primdecZ (cf. 
section|3]) and its parallelization implemented in Singular. Timings are conducted 
by using the 32-bit version of Singular 3-1-2 on an AMD Opteron 6174 with 48 
CPUs, 800 MHz each, 128 GB RAM under the Gentoo Linux operating system. 
All examples are chosen from The SymbolicData Project (cf. [G]). 

Remark 4.1. The parallelization of our algorithm is attained via multiple processes 
organized by Singular library code. Consequently a future aim is to enable par- 
allelization in the kernel via multiple threads. 

Remark 4.2. In Singular one can compute Grobner bases not only over fields 
but also over the rings Z and Z/mZ (resp. Z/2'Z as a special case of Z/mZ). 
For the integers the implementation is based on the theory for Grobner bases over 
integral domains as introduced by Adams and Loustaunau (cf. |AL) . chapter 4). 
For factor rings further theory needed to be developed by the Singular- Team in 
Kaiserslautern (cf. |GSW] . |Wj ). Details about the corresponding implementation 
are presented by Wienand (cf. [W]. chapter 3). 
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We choose the following examples: 

Example 4.3. Coefficients: integer, ordering: drf^j Gerdt-93a.xml (cf. [G ) con- 
sidered with another integer generator 2-3-5-13-17-181. 

Example 4.4. Coefficients: integer, ordering: dp, Gerdt-93a.xml (cf. [TJ]) con- 
sidered with another integer generator 2-3-5-13-17-31-181. 

Example 4.5. Coefficients: integer, ordering: dp, Gerdt-93a.xml (cf. [G]) con- 
sidered with another integer generator 2 ■ 3 • 5 • 13 - 17 - 31 - 37 - 181. 

Example 4.6. Coefficients: integer, ordering: dp, Steidel_6 . xml (cf. ES ) con- 
sidered with another integer generator 2 • 3 • 5 • 7 • 11 • 13 • 17 • 19 • 23. 

Example 4.7. Coefficients: integer, ordering: dp, Steidel_6 . xml (cf. ES ) con- 
sidered with another integer generator 2 • 3 2 • 5 • 7 3 • 11 • 13 • 17 • 19 • 23. 

Example 4.8. Coefficients: integer, ordering: dp, Gonnet-83 . xml (cf. |BGKj ) 
considered with another integer generator 2-3-5-7-11-13-17-19-23. 

Example 4.9. Coefficients: integer, ordering: dp, Gonnet-83 . xml (cf. [BGK ) 
considered with another integer generator 2 2 • 3 2 ■ 5 • 7 ■ 11 • 13 • 17 ■ 19 ■ 23. 

Table Q] summarizes the results where primdecZ* (fc) denotes the parallelized 
version of the algorithm using k processes. All timings are given in seconds. 



Example 


primdecZ 


primdecZ* (2) 


primdecZ* (3) 


primdecZ* (4) 


H3I 


604 


383 


339 


249 


14.41 


757 


480 


392 


350 


14.51 


907 


542 


396 


396 


14.61 


17 


9 


7 


4 


14.71 


10 


6 


5 


4 


14.81 


21 


14 


10 


8 


14.91 


39 


35 


34 


31 



Table 1 . Total running times for computing a primary decompo- 
sition of the considered examples via primdecZ and its parallelized 
variant primdecZ*(fc) for k = 2,3,4. 
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Degree reverse lexicographical ordering: Let x a ,x@ be two monomials in x, i.e. a, /3 E 
N n . x a > dp xP deg(a; a ) > deg(xP) or (deg(x a ) = deg(xP) and 31 < i < n : a n = 

/?„,..., oci-i = /8i_i,aj < Pi), where deg{x a ) = ai + . . . + a n ; cf. |GP| . 
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